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Abstract: We re-evaluate the prompt atmospheric neutrino flux, using the measured 
charm cross sections at RHIC and the Large Hadron Collider to constrain perturbative 
QCD parameters such as the factorization and renormalization scales, as well as modern 
parton distribution functions and recent estimates of the cosmic-ray spectra. We find 
that our result for the prompt neutrino flux is lower than previous perturbative QCD 
estimates and, consequently, alters the signal-to-background statistics of the recent IceCube 
measurements at high energies. 
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1 Introduction 

Decays of mesons produced when cosmic ray primaries at high energies interact with nuclei 
in the atmosphere lead to the production of a large flux of neutrinos. The resulting neutrino 
flux at energies from 10 GeV-1 TeV has been observed by several pioneering experiments 
[1-8] over the last two decades. The dominant contribution to the neutrino flux at these 
energies comes from the decays of charged pions and from leptonic 

and semi-leptonic decays of kaons [9]. Of late, the construction of first the AMANDA 
detector [10] and then, its upgrade to the first > 1 km 3 volume neutrino telescope, the 
IceCube (IC) Neutrino Observatory, has made the detection of neutrinos at much higher 
energies E u > 100 TeV possible (see [11] for a review). Moreover, there are plans for the 
KM3NeT detector in the Mediterranean [12], and for an IceCube-Gen2 upgrade [13], which 
will reach even higher neutrino energies. While AMANDA had already seen atmospheric 
neutrinos in the high TeV’s [14], IC has recently revealed its first observation of high energy 
neutrinos with energies going up to 2.1 PeV [15, 16]. In total, IC observes 37 events in 
the energy range 10 TeV-2.1 PeV over 988 days, with atmospheric neutrinos and muons 
expected to be responsible for about 10-11 of those, their fluxes being especially significant 
for energies up to 200 TeV. A precise understanding of the prompt neutrino background 
at these energies is crucial in determining the IC signal statistics. Because present values 
for the flux used by IC depend on theoretical results 1 [18-22], it is timely to revisit the 
theoretical predictions in light of updated information on the various QCD inputs as well 

TceCube has recently published an upper limit [17] on the prompt flux at energies above 1 TeV at less 
than 1.52 times the flux predicted in ref. [18]. 


- 1 - 





as recent high-energy experimental results such as those obtained at the Large Hadron 
Collider (LHC). 

Prompt neutrinos, produced in the atmosphere by decays of charmed hadrons ( e.g., 
D°, D ± —> u e ^X) that come from cosmic-ray interactions with atmospheric nuclei (pN —> 
ccX ), become an important component of the atmospheric neutrino flux at these high 
energies [18-22], For the IC study of neutrinos from astrophysical sources, atmospheric 
neutrinos are a background and therefore, need to be carefully estimated. The neutrino flux 
from pions and kaons, hereafter called the conventional neutrino flux, has been previously 
calibrated against multi-GeV atmospheric neutrino observations. This conventional flux is 
well understood, but the prompt neutrino flux, being sub-dominant at lower energies, is 
less so. 

The evaluation of the prompt flux requires as a first step the computation of the charm 
production cross-section in the process pN —> ccX, followed by the hadronization of charm 
particles. The computation of charm production in hadronic collisions introduces some 
of the biggest uncertainties that translate into the evaluation of the prompt flux. These 
arise due to uncertainties in the knowledge of a) the charm mass ( m c ), b ) the factorization 
(Mp) and renormalization (pr) scales, and c) the choice of parton distribution functions 
(PDF’s). In addition, the final step in computing the prompt lepton flux involves folding 
the charmed hadron production cross-section with the flux of incoming cosmic-rays, and 
thus incurs large uncertainties from the limited understanding of the extremely high-energy 
cosmic ray composition. 

Two approaches to the calculation of the prompt lepton flux are via perturbative QCD 
(pQCD) in the parton model (as in, e.g., [20-22]) and with the dipole model (see, e.g., 
[18]). Recently, the production of charm mesons in the atmosphere has also been studied 
using the event generator SIBYLL [23]. In this paper we focus on a new evaluation of the 
prompt flux using next-to-leading order (NLO) QCD. We use recently updated PDF’s in our 
pQCD calculations while limiting the uncertainties in the charm production computation 
by incorporating constraints from related high energy collider results. Finally, we evaluate 
the prompt neutrino flux and, therefore, the prompt atmospheric neutrino background for 
IC by using recent estimates of the cosmic-ray spectrum and composition [24-26] and our 
updated charm production cross-section. 

Our focus here is on the link between the charm cross section measurements, including 
LHC results, and the prompt neutrino flux in the context of pQCD. This allows us to 
make estimates of some of the uncertainties in the flux evaluation. This work is part of 
a broader program to assess the uncertainties in the theoretical evaluation of the prompt 
flux accounting for a wider range of approaches to calculating the charmed hadron cross 
sections [27]. 

In the next section, we review the pQCD evaluation of the production of charmed 
hadrons. In section 3, we describe the Z-moment approach to evaluation the flux, and we 
show our results for the prompt muon (and electron) neutrino plus antineutrino fluxes. We 
briefly discuss the implications of this new flux for the background evaluation at IceCube 
in section 4, before drawing our conclusions in the final section. 
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Expt. 

\fs [TeV] 

a [mb] 

PHENIX [31] 

0.20 

0.5611°;*? (sys) 

STAR [32] 

0.20 

0.797 ± 0.210 (stat)l°;*« (sys) 

ALICE [28] 

2.76 

4.8 ± 0.8 (stat)!^ (sys) ± 0.06 (BR) 
± 0 . 1 (frag) ± 0.1 (lum)lg f [ (extrap) 

ALICE [28] 

7.00 

8.5 ± 0.5 (staQlg 4 (sys) ± 0-1 (BR) 
±0.2(frag) ± 0.3 (lum)!^ (extrap) 

ATLAS [29] 

7.00 

7.13 ±0.28 (stat)18;88 ( s y s ) 

±0.78 (lum)l^ ( e xtrap) 

LHCb [30] 

7.00 

6.100 ±0.930 


Table 1: Total cross-section for pp(pN) — > ccX in hadronic collisions, extrapolated based 
on NLO QCD by the experimental collaborations from charmed hadron production mea¬ 
surements in a limited phase space region. 

2 Charm production cross section 

The PeV energy range for atmospheric neutrinos corresponds to an incident energy E p ~ 
30 PeV for pA fixed target interactions. The LHC center of mass energy y/s = 7 TeV 
is equivalent to a fixed target beam energy in pp collisions of Eb = 26 PeV. The LHC 
measurements of the charm production cross section [28-30] together with recent RHIC 
[31, 32] and modern parton distribution functions (PDFs) have narrowed down some of the 
uncertainty in the rate of charm production in the atmosphere. The experimental results 
at high energy for the charm production cross-section in hadronic collisions are listed in 
Table 1. 

In Ref. [33], Nelson, Vogt and Frawley have investigated a range of factorization and 
renormalization scales using the CT10 PDF’s [34] and the NLO order QCD calculation of 
Nason, Dawson and Ellis [35, 36]. Using a charm quark mass central value of m c = 1-27 
GeV based on lattice QCD determinations of the charm quark mass, as summarized in 
Ref. [37], and a combination of fixed target, PHENIX, and STAR charm production cross- 
sections, they find that AIp/m c = 1.3-4.3 and PR/m. c = 1.7-1.5 with Mp = 2.1 m c and 
Pr = 1.6m c as central values. We use these values of parameters as a guide to the range 
of theoretical NLO charm cross sections expected at high energies. 

In our calculation we use the NLO Fortran code of Cacciari et al. [38, 39] that includes 
the total cross section [35] as well as the single [36] and double differential [40] distributions 
of charm (i.e., da/dy and d 2 cr/dydpT respectively). The cross sections shown in figure 1 for 
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proton interactions with the iso-scalar nucleon 2 pN —> ccX are obtained using factorization 
and renormalization scales relative to the charmed quark transverse mass rn^ = (m 2 + Pp) 
rather than relative to the charm quark mass m c . This choice of the scale dependence is an 
acceptable representation of the factorization and renormalization scales both at moderate 
energies (where the choice of scale is not sensitive to pp) and at high energies (where the 
Pt dependence in the scales becomes important). The CT10 PDFs are our default choice 
for the NLO evaluations of charm production. The shaded band in figure 1 comes from 
using (Mp, pr) = (1.25, 1.48) mp and {Mp, pp) = (4.65, 1.71) nip (labeled as CT10 
Limits), as suggested in ref. [33]. The central solid curve (CT10 Central) is obtained using 
{Mp, pp) = (2.1, 1.6) mp. The LHCb [30], ALICE [28] and ATLAS [29] cross sections 
measurements are within our calculated cross section uncertainty band, thus we use this 
range of scales for our flux uncertainty band. The range of QCD parameters consistent with 
observations have been studied in [33] by taking into consideration both total as well as 
differential cross-section measurements at the LHC for a wide range of the kinematic phase 
space. This includes measurements at low transverse momenta and forward rapidities. The 
upper and lower limits of our QCD calculation, consistent to within the error-bars for the 
experimental data, are based on these parameter uncertainties (see e.g., shaded area in 
figure 1). 

The total pp —> ccX cross-section at high energies is dominated by gluon fusion, so 
the cross section is sensitive to the gluon PDF small-a: behavior. The perturbative results 
in [21] overestimate the charm cross-section because of the steeply rising gluon PDFs at 
small-® for the CTEQ3 PDF sets [41] used, e.g., in [18, 21], By comparison, the updated 
CT10 PDFs we use here have a more slowly growing behavior as x decreases. These PDFs 
are fit based on recent experimental inputs (including from LHC, Tevatron and others). 
Using the MSTW [42] PDFs, which have an even smaller gluon PDF growth at small-®, 
leads to cross-sections that grow even more slowly (compared to CT10) at higher energies, 
however, the results still fall within the theoretical uncertainty band shown in figure 1. 

By choosing Mp and pp with mp —> m c , we find a similar energy dependence of 
the total charm cross section [33]. In addition, the central values for the charm cross 
sections are roughly comparable for the same choice of PDF. Using the CT10 PDF and 
m c = 1.27 GeV, in the energy range Ej, = 10 6 — 10 8 GeV, the cross section obtained 
with (Mp , pr) = (2.1, 1.6) m c is within about ±10% of the cross section calculated with 
{Mp, pp) = (2.1, 1.6) mp. The upper choice of scale factors multiplying m c gives a 
comparable cross section for the same factors multiplying mp. For the lower range of scale 
factors, the high energy charm pair production cross section is smaller for m c dependent 
scales. For example, for E}, = 10 s GeV using ( Mp , pr) = (1.25, 1.48) m Cl the cross section 
is about 40% lower than the cross section evaluated using mp dependent scales as shown 
by the grey fine-dashed curves in figure 1. Consequently, while the central values for either 
choice of the scale give almost identical prompt neutrino flux, it is the uncertainty band 
that is smaller for the ~ mp scale choice. 

As discussed below, we use the superposition approximation in describing the cosmic 

2 Since the charm production cross section is dominated by gluon fusion, p ~ n ~ N. 


- 4 - 




Figure 1: The charm production cross section a p N->-cc+x at NLO with rn c = 1.27 GeV 
using the CT10 parton distributions for a range of scales described in the text, with the 
central set with factorization and renormalization scales Mp = 2.10 mp and p,p = 1.6 
respectively. Apart from experimental data points listed in table 1, results from HERA-B 
[43] and lower energy experiments summarized in [44] for pN scattering are shown (labelled 
as Fixed target expts .). For comparison, we also show the lower and upper limits (grey 
fine-dashed curves) when the renormalization and factorization scales are made to vary 
proportionally to rn c rather than to mp. 


ray interactions with air nuclei [45]. In this approximation, the cosmic ray all-particle flux 
is reduced to a nucleon flux. We require the NA —>• ccX cross section, and use p ~ n ~ N 
(as noted previously). We continue to label the cosmic ray nucleon with p, to distinguish it 
as the beam nucleon as opposed to the target nucleon. The average atomic number of the 
target air nuclei is A = 14.5. Measurements of charm production on a variety of targets 
[46, 47] show that a{jpA —> ccX ) ~ Aa(pN —> ccX). We use this approximation here. 

2.1 Differential cross section 

While we seek compatibility of the total charm quark pair production cross section with 
the results reported by the experimental collaborations, the dominant contribution to the 
prompt flux is from forward production of charm, including fragmentation into charmed 
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hadrons. In our semi-analytic evaluation of the prompt atmospheric lepton flux, we require 
the differential charmed hadron energy distribution, 


da 

dE h 



kX)D h k 



dE k 

E k 


( 2 . 1 ) 


in terms of the parton level differential distribution and the fragmentation function D k . 
Here, h = D , D^r. and k = c, c. We approximate the fragmentation functions 

for charmed hadrons as energy independent. Eq. (2.1) can be written as 

~T—{pA hX) =A f 1 (pN -> cX)D h c {z) (2.2) 

dx E J XE z dx c 

in terms of xe = Eh/Ei, and x c = E c /Eb = xe/ z for an incident cosmic ray nucleon energy 
(beam energy) E In figure 2, we show the differential cross-section as a function of x c 
for Eb = 10 3 , 10 6 and 10 9 GeV in pN scattering, here for (Me, Hr) = (2.1, 1.6) my- The 
distributions for (Mr, Pr) = (2.1, 1.6) m c are very similar. 

We can compare our results here to those obtained previously, notably in [21], With 
the CT10 NLO PDF’s, the uit dependent scales, and a full NLO calculation, we find that 
our differential distribution at low x is lower than in ref. [21] at high energies ( e.g ., about 
28% lower at 10 9 GeV for x = 0.1). As previously discussed, this stems from the use of 
updated PDFs which have a slower growth at small x than the CTEQ3 PDFs used in [21]. 

Our default choice of fragmentation functions for charmed hadrons is that of Kniehl 
and Kramer [48]. The net effect of including the fragmentation functions is to reduce the 
predicted flux by about ~ 30% relative to the flux without fragmentation included. 


3 Prompt lepton flux 


We use the Z-moment approach [49, 50], including an energy dependence of the Z-moments 
[20] and approximating the depth of the atmosphere as infinite. In the exponential atmo¬ 
sphere approximation where the density is 


p(h) = po exp(—/i//ro) 


(3.1) 


for po = 2.03 x ICG 3 g/cm 3 and ho = 6.4 km, the low energy and high energy lepton 
fluxes have particularly simple forms, involving the spectrum weighted Z-moments. The 
production moments are defined by 

rl dXR 4>p(E h /xE) 


Zph(E h ) = 


1 


.da . Ar 
x A- — (pN 


' xe 


x E 4>p(E h ) a pA {E h ) dx E 


hX) . 


(3.2) 


The all-nucleon cosmic ray flux as a function of atmospheric column depth X is 


<t>p(E, X) ~ 4>° p {E) exp(-V/A p ) = (dN/dE) exp(- X/A p ) , (3.3) 

with A p = X p (E)/(l — Zpp(E)). For the proton-air cross section, we use an approximate 
parametrization of the EPOS 1.99 cross section [51] that is consistent with the high energy 
results of the Pierre Auger Observatory [52], 
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Figure 2: The differential cross section da/dx c for the charmed quark, as a function of x c 
for E = 10 3 , 10 6 , 10 9 GeV for m c = 1.27 GeV and Mp = 2.1 my, fiR = 1.6 mp using the 
CT10 NLO PDFs. 


For decays, the differential cross section is replaced by the differential decay distribu¬ 
tion, and the cosmic ray flux is replaced by the charmed hadron flux. At high energies, we 
evaluate the high energy decay Z-moment with a spectral weight of (f)p(E / x e) / <j>p{E) since 
4>h oc 4> p at high energies. The low energy decay Z-moment is evaluated using a spectral 
weight of 1 /xe- since the low energy charmed hadron flux is proportional 

to EcjP p {E). 

For low energy, the lepton l = /i, v t flux from h —> £ decays is approximated by 

( 3 . 4 ) 

App 


while for high energy (see e.g. [20]), 


$ i9h (h) = Z ( 


_ yhigh Zph ln(Afc/ A p ) q 


hi 


1 Z p p 1 Ap/A/j 


(3.5) 


Each Z factor, A and has an energy dependence, suppressed in the notation above. The 
resulting lepton flux from charmed hadrons h is 


h 


<P l e ow (h)cj> h e i9h (h) 
<$ w (h) + <l>f lh (h) ' 


(3.6) 


The lepton fluxes from atmospheric charm depend on the cosmic ray flux directly 
through 4>p(E) and in the evaluation of the energy dependent Z-moments. 
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3.1 Cosmic ray flux 


The cosmic ray flux has been measured directly and indirectly over a wide energy range. 
Direct measurements are available to energies of about ~ 100 TeV. At higher energies, 
indirect measurements are made by air shower array experiments. While there are some 
discrepancies in the normalization of the cosmic ray spectrum at high energies, overall the 
all-particle cosmic ray spectrum for the energy range of interest, 10 3 -10 10 GeV, approxi¬ 
mately follows a broken power-law with the break occurring at E ~ 5 x 10 6 GeV. Many 
earlier evaluations of the prompt lepton flux [18, 20-22, 53] used the broken power-law 
form for the nucleon flux with [ 20 ] : 


$(E) 


1.7 E~ 2 - 7 for E < 5 • 10 6 GeV 
174 E~ 3 for E > 5 • 10 6 GeV, 


(3.7) 


for E in GeV and the nucleon flux in units of cnr 2 s _1 sr -1 GeV ^ 1 . With the fairly recent 
measurements from ATIC [54], CREAM [55, 56] and Pamela [57], combined with earlier 
measurements, Gaisser [24] and collaborators [25, 26] have taken a multicomponent model 
with three or four source populations to develop models for the cosmic ray composition. 
Their parametrizations depend on the particles’ electric charges Z and maximum energies 
of the source populations, with spectral indices 7 that vary by population and nucleus. We 
use here the parametrization by Gaisser in ref. [24] with three populations: from supernova 
remnants, from other galactic sources and from extragalactic sources. The H3a flux from 
ref. [24] has a mixed composition in the extragalactic population, while the extragalactic 
population in what we call the H3p flux is all protons. Thus, the cosmic ray nucleon 
spectrum is identical for H3a and H3p for nucleon energies below ~ 10' GeV. Converting 
the all-particle flux to the nucleon flux, the H3a and H3p fluxes are shown along with the 
broken power-law in figure 3. 

The composition of the cosmic rays causes a much steeper drop in the nucleon flux 
above the knee energy than when using the simple broken power-law parametrization. This 
is particularly important for the high energy prompt lepton flux. To allow for comparisons 
with earlier work, we show our results for the prompt lepton flux for the broken power-law 
and the H3a and H3p cosmic ray fluxes. 


3.2 Z-moment and prompt lepton flux results 

The production Z-moments are shown as a function of energy in figure 4a. For each of Z pD 0 
and Z pD ±, the three curves show the moments evaluated for the three respective cosmic-ray 
nucleon fluxes presented in figure 3. For the H3p flux, we also show the band of Z-moments 
from the range of differential cross sections by taking (Mp, fip) = (1.25, 1.48) mp (for 
lower limit) and (Mp, fip) = (4.65, 1.71 )mp (for upper limit). This relative band of 
variation is identical for the other Z-moments shown in the figure. Figure 4b shows the 
ratio of the central Z p£) o-moments obtained using the H3a and H3p fluxes to that evaluated 
using the broken power-law nucleon flux. 

The major difference between the D-meson production Z-moments when using the 
power-law CR flux [from eq. (3.7)] against a more recent CR flux estimate, such as the 
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Figure 3: The all-nucleon cosmic ray spectrum as a function of energy per nucleon for 
the three component model of ref. [24] with a mixed extragalactic population (H3a) and 
all proton extragalactic population (H3p), and for the broken power-law of eq. (3.7). 


Gaisser H3p flux, arises at the high energies ^ 10' 5 GeV, where the latter noticeably dip, 
before rising sharply at energies beyond the tens of PeV. In contrast, the Z-moments when 
using the broken power-law follows a more steady behavior. This difference in nature can 
be traced to the particular behavior of the Gaisser cosmic ray primary fluxes—a significant 
softening of the spectral shape occurs at around a few PeV energies, where the population 
transitions from being dominantly galactic to extra-galactic, before the spectra hardens 
again at energies around a few hundred PeV (see figure 3). When translated to the pro¬ 
duction Z-moments, these effects are visible at comparatively lower energies because of the 
inelasticity of the high energy pp collision, which implies that only a small fraction (given 
by (xe) ~ 0.1) of the incident proton energy goes into the produced cc. The nature of 
the Z-moments, in turn, translates directly to the total prompt lepton flux (as shown in 
figure 5a). The central Z-moments obtained using the H3p estimate will henceforth be our 
benchmark result when determining the prompt flux and correspondingly the event-rates 
at IC. 

As discussed above, we use the charmed hadron spectral weights for the decay Z-moments. 
These are evaluated using dn/dE from ref. [49, 50, 58]. 

Additional Z-moments are needed for the flux evaluation, in particular Z pp and Zhh 
along with A h- For Z pp , we have approximated the pA —» pX differential cross section with 
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(a) (b) 

Figure 4: (a) Production Z-moments for pA M for M = D° + D° (as thick curves) 
and M = (thin curves) for H3a (green dot-dashed curves), H3p (orange solid curves) 
and broken power-law (blue short-dashed curves) cosmic ray fluxes. The range of variation 
for the Z-moments is shown (orange shaded region) for the Z pE o when the H3p cosmic ray 
flux is used for computation, and the relative range of variation is identical for the other 
curves. 

(b) Ratio of central Z-moments for pA —> D° + D° using the Gaisser H3a and H3p fluxes 
to the broken power-law cosmic ray flux (from eq. (3.7)). 

a scaling form 

^~a pA (E)(l + n)(l-x E ) n (3.8) 

with <JpA as described above and n = 0.51. With these choices, at E = 10 3 GeV for the 
broken power-law, Z pp = 0.271 and A p = A p /(1 — Z pp ) = 116 g/cm 2 . By comparison, 
the scaling values in [49, 50] are Z pp = 0.263 and A p = 117 g/cm 2 . The change in the 
cosmic ray spectrum with the broken power-law together with the energy dependence of 
the cross section reduces our calculated Z pp to 0.231 and A p to 67 g/cm 2 at E = 10 s GeV. 
Similar results are obtained for the H3a and H3p cosmic ray flux inputs. We remark that 
in [20], energy dependent Z-moments evaluated using the PYTHIA Monte Carlo program 
[59] are used, giving e.g. Z pp ( 10 3 GeV) ~ 0.5. These are also used in [18]. The low energy 
flux is proportional to (1 — Zpp)^ 1 , so this numerical factor is important to the overall 
normalization. 

For the charmed hadrons’ interaction lengths and interaction Z-moments, we use kaon- 
proton interactions as representative. For all charmed hadrons, we use the same expres¬ 
sions, based on kaons. We take 

- A°- 75 a KN (E)( 1 + n)(l - x E ) n (3.9) 

ax e 

with (Ten determined by the COMPAS group and summarized by the Particle Data Group 
in [37]. We find that setting n = 1 gives Zkk = 0.217 and An = 162 g/cm 2 for the broken 
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power-law at 10 3 GeV, reducing to Zkk = 0.176 and A k = 40 g/cm 2 at 10 8 GeV. The 
scaling values in [49, 50] are 0.211 and 175 g/cm 2 , respectively. The precise value of n for 
meson scattering in eq. (3.9) affects only (fy l9h . 

Our prompt lepton fluxes are shown in figure 5a. We show E s (j) u +& as a function of 
neutrino energy. The fluxes of // + p, and v e + v e are the same as shown in the figure, in 
the approximations we use here. The upper band shows our NLO result using CT10 PDFs 
with the range of (Mp, hr) discussed in section 2, using the broken power-law as the input 
cosmic ray all-nucleon spectrum. The lower band shows the prompt flux using the H3p 
cosmic ray flux inputs. The H3a cosmic ray flux results in a lower prompt lepton flux for 
energies above ~2x 10 5 GeV, roughly a factor of two lower at E = 10 8 GeV. 




(a) (b) 

Figure 5: (a) Our benchmark results for the prompt flux scaled by E 3 is shown as 

an orange curve, with the cosmic ray flux given by the Gaisser H3p fit (see figure 3). The 
blue curve uses instead a broken power-law (as used in previous analyses, e.g., [18]). For 
each curve, the associated shaded region indicates the uncertainty due to variation of the 
QCD parameters. The vertical conventional flux from Honda (see, e.g., [20]), reweighted 
to the H3a cosmic-ray primary flux, is also shown. 

(b) Comparison of neutrino fluxes with variation in scales and PDFs for the broken power- 
law CR primary flux. Shown are the results for central values obtained using CT10 as the 
PDF with (Mf,/Ur) oc mx (solid slate blue line) and with scales (Mp, hr ) oc m c (dashed 
slate blue curve) along with their associated bands of variation (corresponding to QCD 
parameters discussed in text) as solid and hatched fills for scales proportional to wt and 
m c respectively. The central flux (corresponding to Mp = 2.1 m c , HR = 1.6?tit) evaluated 
for CTEQ3 as the PDF is shown as the pink dashed curve, along with the dipole model 
computation (gray short-dashed curve) of ref. [18] . The flux uncertainty from [18] is shown 
as a grey band. For comparison, the vertical conventional flux from Honda (see, e.g., [20]), 
based on the broken power-law cosmic-ray primary flux, is also shown. 
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3.3 Comparison with previous prompt neutrino flux calculations 

In [18], the dipole model approach including effects of parton saturation was used for charm 
production in pp —> ccX. With respect to the total neutrino flux obtained therein, we And 
that our benchmark results, i.e., obtained using Gaisser H3p cosmic ray flux as opposed 
to the use of broken power-law spectra in the former, yields prompt neutrino fluxes that 
are reduced by a factor which varies from about 2 at lower energies (below 100 TeV) to 
a maximum reduction factor of about 8 at high energies (few PeV). On the contrary, the 
result for the central value of the flux obtained here using a broken power-law cosmic ray 
spectra, for comparison, is seen to be lower compared to the results in [18] by up to a 
factor of 3. This reduction, despite the use of the same cosmic ray spectra, mainly arises 
due to the updated evaluation of the Z pp moment [as explained just below eq. (3.8)], and 
from differences in the large xe values of da/dxE in the dipole and perturbative QCD 
approaches. Figure 5b additionally shows the flux that arises for the central scale choice 
performed with the CTEQ3 PDFs and the broken power-law in eq. (3.7). We And that the 
steep increase in cross-section with energy that one obtains when using this older PDF set, 
translates directly to larger Auxes at energies of 10 TeV and higher (compared to central 
CT10 fluxes obtained using the same cosmic-ray primary flux); nonetheless, it turns out to 
be smaller at low energies than that obtained in the dipole computation in [18]. Finally, 
figure 5b also shows that using scales proportional to m c rather than m t leads to similar 
central predictions but somewhat larger bands of uncertainty. This is a feature that can 
be traced back to the pp —> ccX cross-section, as discussed in section 2. 

For both figures 5a and 5b, we show the conventional atmospheric neutrino flux from 
ref. [60] . Since the original computation of this flux in [60] was based on the broken power- 
law cosmic ray primary flux, for consistency, we reweigh this flux in figure 5a to account 
for the updated H3a cosmic-ray primary flux 3 using the code NeutrinoFlux [61]. This 
reweighted flux is used also to determine the resulting conventional atmospheric background 
in figure 6. 

Quantifying, approximately, the changes in our evaluation over that in [18] aside from 
using perturbative QCD rather than the dipole model approach, we note that the following 
factors are responsible for a reduced prompt flux: 

• Use of recent cosmic-ray flux estimates which predict a lower nucleon flux at high 
energies (10 7 -10 9 GeV) than the previous broken power-law (see eq. (3.7)). Cor¬ 
respondingly, the prompt flux evaluated using the H3p cosmic ray primary flux is 
reduced by a factor of 1.2-3 (figure 5a) relative to the broken power-law between the 
energies 10 5 -10 6 GeV, for example, where the prompt background is expected to be 
important. 

• Reevaluation of Z pp and A p , and Z^ and A x , using the H3p cosmic ray flux and 
a parametrization of the p —> p, D —> D ~ K —>• K differential cross sections, and 
using an updated p-Air cross-section, decreases the flux by about 30%. 

li The difference between results obtained using H3a and those using H3p fluxes is negligible below 1 
PeV. Thus, we do not distinguish between the two when computing the conventional flux, which is only 
significant at much lower energies, i.e., up to 200 TeV. 
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• Use of CT10 PDFs instead of the older CTEQ3 PDFs used in [21] reduces the flux 
by an additional factor of 2-3 at the high energies (> 10 4 5 GeV), as figure 5b shows 
by comparing fluxes computed using the same broken power-law cosmic-ray primary 
flux and QCD parameters, but with the two different PDF sets. 

• Between the dipole calculation in [18] and the results shown here for the broken 
power law cosmic ray flux, the large x behavior of the dipole and perturbative QCD 
evaluations of pA —> ccX, input to the evaluation of the production Z-moments, 
account for the relative fluxes differing by a factor of ~ 1.5. Uncertainties associated 
with variations of the QCD parameters for the perturbative and dipole evaluations 
are shown with the (overlapping) bands in figure 5b. 

4 Prompt background at IceCube 

The IceCube atmospheric neutrino background event rates, both prompt and conventional, 
can be evaluated from the corresponding flux by making use of the flavor-dependent ef¬ 
fective areas for the detector, given in [15]. The central value of the estimated prompt 
neutrino background relevant to the 988-day IC observation is shown in figure 6. Addi¬ 
tionally, the figure also shows the maximum and minimum estimates due to variation in 
the cross-section (see figure 1) and differential charm quark energy distribution obtained 
by varying the QCD parameters compatible with experimental results. 4 We note that the 
procedure using the IC effective areas overestimates the prompt event rates because it does 
not take into account the neutrino self-veto as discussed in [62]. 

Using the Honda estimates for conventional atmospheric neutrino flux, 5 events from 
background muons estimated in [16], and our central prompt flux, we compute the total 
atmospheric background. This is shown in figure 6 as the blue dashed curve, with the 
uncertainty indicated as a hatched area around this central curve. When compared with the 
conventional atmospheric neutrino background, the prompt flux is only a minor contributor 
at energies where the background is significant (about 200 TeV and lower). Consequently, 
the total atmospheric neutrino background at these energies is not significantly affected by 
the ~ 10%-15% uncertainty in the prompt flux from variation of the QCD parameters. 

It is clear from a comparison with the total atmospheric background (gray dashed 
curve) in figure 6, obtained using a previous results for the prompt neutrino flux [16] 
and IC’s experimental 90% CL prompt flux limits, 6 that our new pQCD results lead to a 
significant reduction in the total background in the lower energies up to 100 TeV. This will 
in turn lead to revisions in the evaluation of statistical significances of the astrophysical 

4 When comparing event rates from prompt neutrino flux obtained here to those obtained in [16], the 
difference is not as stark as one would expect from a straight comparison at the level of fluxes obtained 
here to that obtained in [18], because IC implements a reweighing procedure on the central estimate of the 
prompt neutrino flux in [18] to already account for updated cosmic-ray spectra. 

5 For the conventional atmospheric neutrino flux too, IC uses the Honda flux reweighted to account for 
changes in fits to the cosmic-ray spectra from older power-law estimates to more recent estimates, e.g. in 
[24]. We have used the same method in figure 5a to obtain the conventional neutrino flux. 

6 IC’s 90% CL limits on the prompt flux corresponds to 3.8 times the central prompt flux estimated in 
[18]- 
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Figure 6: Event rates at IceCube from prompt neutrinos, with our updated prediction for 
the prompt flux indicated in magenta, along with uncertainties from variation in the QCD 
parameters indicated as a hatched region around the central curves. The re-evaluated to¬ 
tal atmospheric neutrino background (blue dashed curve) includes events from our central 
prompt prediction and the Honda conventional neutrino flux reweighted for H3a cosmic ray 
flux (with the uncertainty in prompt event rates indicated as a hatched region around it). 
The total neutrino background estimated by IC at the level of 90% CL charm limit [16] is 
shown (dashed gray curve) for comparison. The prediction for total event rates using the 
E~ 2 fit astrophysical signal from [16] and updated atmospheric background (i.e., includ¬ 
ing our re-evaluated prompt background) is shown as a green curve, while the similar IC 
estimate using the older atmospheric background and prompt at the level of 90% CL (see 
[16]) is shown as the gray thick-dashed curve. For both these latter curves, the background 
includes contribution from atmospheric muons (reproduced from [16]), in addition to at¬ 
mospheric neutrinos. Observed total event rates at the IC are shown as solid red blocks, 
along with their associated lcr statistical uncertainty. 


signal for the IC, but a detailed re-analysis of the IC signal is beyond the scope of the 
present work. We note that this reduced background implies that the IC fixed-slope E~ 2 
best-fit astrophysical flux underestimates the low-energy signal in the observed events. 
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5 Conclusions 


Hadronic interactions at extremely high energies necessarily involves understanding of QCD 
at these energies; however, this is beset by relatively large uncertainties in the related 
computation. Interactions between cosmic ray protons and atmospheric nuclei which lead 
to the production of neutrinos via the formation of a mesons (tt^, Lf 0, , D 0,± , D±) 
and baryons, and its consequent decay is one such process. In the era of the IceCube 
observations of ultra-high energy neutrinos, the importance of properly understanding and 
estimating the background from atmospheric neutrinos cannot be overstated. Since earlier 
perturbative QCD results for prompt neutrino fluxes were based on fairly old PDFs and 
were not constrained by the recent LHC data, we have revisited the computation in the 
present work, incorporating recent developments in the understanding of low-x PDF’s and 
using recent experimental results that have bearing on the relevant QCD inputs. 

By varying the QCD parameter space including the involved factorization and renor¬ 
malization scales within bounds set by results from LHCb, ALICE, ATLAS among others, 
and making use of recent CT10 PDF’s, we have re-evaluated cr(pp —> ccX). We have de¬ 
termined a(pN —> ccX ) from the latter, and thereafter, using the standard procedure of 
energy-dependent Z-moments, we have computed the total prompt flux expected at the 
detector. For an estimate of the composition and spectrum of the cosmic ray primaries, we 
have relied on the recent results from ref. [24] rather than use a broken power-law which is 
now known to be an inaccurate representation of the cosmic ray nucleon flux at the high 
energies. 

We find that the prompt neutrino background at the IceCube is lower than that esti¬ 
mated in [15]. We note that this may lead to a consequent revision in the estimates of the 
statistical significance of the IC signal over the now reduced atmospheric background. 
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